PROBABILITY AND FOURIER DUALITY FOR AFFINE ITERATED 

FUNCTION SYSTEMS 



DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 

Abstract. Let d be a positive integer, and let /x be a finite measure on M. d . In this paper 
we ask when it is possible to find a subset A in R d such that the corresponding complex 
exponential functions e\ indexed by A are orthogonal and total in £ 2 (/i). If this happens, 
we say that (/x, A) is a spectral pair. This is a Fourier duality, and the x- variable for the 
L 2 (/i) -functions is one side in the duality, while the points in A is the other. Stated this 
way, the framework is too wide, and we shall restrict attention to measures fi which come 
with an intrinsic scaling symmetry built in and specified by a finite and prescribed system 
of contractive affine mappings in R d ; an affine iterated function system (IFS). This setting 
allows us to generate candidates for spectral pairs in such a way that the sets on both sides 
of the Fourier duality are generated by suitably chosen affine IFSs. For a given affine setup, 
we spell out the appropriate duality conditions that the two dual IFS-systems must have. 
Our condition is stated in terms of certain complex Hadamard matrices. Our main results 
give two ways of building higher dimensional spectral pairs from combinatorial algebra and 
spectral theory applied to lower dimensional systems. 
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1. Introduction 

The use of traditional Fourier series has up to recently been restricted to the setting of 
Fourier duality between groups; in the abelian case |Rud62j . between compact groups (such 
as tori) on the one side, each group coming with its Haar measure; and discrete abelian 
groups (such as lattices) on the other. However in dynamics and in other applications 
to computational mathematics, one is often faced with sets arising as attractors, highly 
non-linear, and coming equipped with equilibrium measures. This has led to attempts at 
adapting traditional Fourier tools to these non-linear and non-group settings. In this paper 
we address the Fourier duality question for affine iterated function systems. For some of 
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the earlier literature we refer the reader to [JP921 [JP931 lJP9l IJP951 IJP98cl IJP98al IDJ07bl 
ID J07cl lDJ07al ID J07dl IDR071 IDHPS081 iJorM IJKS071 IOS05] . 

Iterated function systems (IFS) in M. d are natural generalizations of more familiar Cantor 
sets on the real line. Like their linear counterparts, they arise as limit sets X for recursively 
defined dynamical systems. While the functions used may be affine, the limit X itself will 
typically be a highly non-linear object, and will include complicated geometries. They arise 
in operator algebras and in representation theory; and they form models for "attractors" in 
dynamical systems arising in nature. For d = 2, the Sierpinski gasket is a notable example, 
and there is a variety of possibilities for d > 2 as well. Each affine IFS X possesses (normal- 
ized) invariant measures /i (see (I2.6P ). naturally associated with the system at hand (denote 
by X the support of the measure). 

Question: When is there some orthogonal basis in the Hilbert space L 2 (X, //) of the form 
{eA | A G A}, where A C M d and e\{x) = exp(2iiix ■ A)? 

Definition. The restricted class of IFSs [i for which a basis {ca}aga can be found are 
called spectral measures, and the functions e\ are said to form a Fourier basis, (/i, A) is then 
called a spectral pair. 

If X is the middle-third Cantor set, Jorgensen and Pedersen }JP98b] proved that there is no 
such Fourier basis. Nonetheless, many spectral measures have since been found within various 
important classes of fractals; and their significance has been explored by many authors. 

Affine IFSs often take the following form: Start with a d by d matrix R and a finite 
subset B in M. d . Then consider the associated set of affine mappings in M. d , Tf, of the form 
Tb(x) = Rr l (x + b), where R is further assumed to be an expanding integer matrix and 
beB cM. d . 

In [DJ07c] we conjectured that the spectral measures arise precisely when there exists 
duality pairing, i.e., another system L such that the two define a complex Hadamard matrix 
H (see equation (11.11) ). the order of the matrix H being the cardinality of B. We proved the 
conjecture in [DJ07cj when an additional assumption, called "reducibility" , is placed on the 
triple (R, B, L). 

Two approaches to IFSs have been popular: one based on a discrete version of the 
more familiar and classical second order Laplace differential operator of potential theory, 
sec [KS W01[ |Kig04[ ILNRC96] : and the second approach is based on Fourier series, see e.g., 



[JP98b, DJ06aJ. The first model is motivated by infinite discrete network of resistors, and 
the harmonic functions are defined by minimizing a global measure of resistance, but this ap- 
proach does not rely on Fourier series. In contrast, the second approach begins with Fourier 
series, and it has its classical origins in lacunary Fourier series |Kah86j . 

Hadamard matrices. Let R be a d x d integer matrix, Bc^ and L C Z d having the 
same cardinality as B, j^B = j^L =: N. We call (R, B, L) a Hadamard triple if the matrix 

7/W >beB,ieL 



is unitary. 

Let fiB be the invariant measure associated to the affine IFS Tf,(x) = R~ x (x + b), b G B. 
(See Theorem 12.41 below) . We conjectured in [DJ07c] that the existence of a set L such that 
(R, B, L) is a Hadamard triple is sufficient to obtain orthonormal bases of exponentials in 
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Conjecture 1.1. Let R be a d x d expansive integer matrix, B a subset of 7L d with G B. 
Let fiB be the invariant measure of the associated IFS (tj^^b- If there exists a subset L of 
Z d such that (R, B, L) is a Hadamard triple and 6 I then Hb is a spectral measure. 

In |DJ07c] we proved that this conjecture is true under a certain "reducibility" assumption. 

Theorem 1.2. Let R be an expanding d x d integer matrix, B a subset of 7L d with G B. 
Assume that there exists a subset L of 7L d with G L such that (R, B, L) is a Hadamard 
triple which satisfies the reducibility condition (see [DJ07c] ). Then the invariant measure hb 
is a spectral measure. In particular the conjecture is true in dimension d = 1. 

In this paper we will show that the conjecture is true also in some other cases, by reducing 
the problem to some "building blocks" in lower dimensions (Theorem I4.4[ Corollary 14. 5p . 
We give an example of an affine IFS when the reducibility condition from |DJ07c] is not 
satisfied but the associated measure is still spectral (Example 15. ip . 



The purpose of the present section is to collect the necessary definitions, and to point 
out the link between the two sides in the general (non-group) Fourier duality; see especially 
Proposition 12.21 and Definition 12.51 Hadamard-triple. 

Definition 2.1. For A G R d , let e x (x) : = e 27riA % x G R d . 

A probability measure on M. d is called spectral if there exists a subset A of lR d such that 
the family of exponential functions {eA | A G A} is an orthonormal basis for L 2 (fi). In this 
case A is called a spectrum for fi. 

The Fourier transform of a measure /i is defined by 



In a number of applications, one encounters a measure jj and a subset A such that the 
functions e\ indexed by A are orthogonal in L 2 (fi), but a separate argument is needed in 
order to show that the family is total. The following is a universal test which applies to any 
subset A: It is a necessary and sufficient condition on a pair (/i, A) allowing us to decide 
whether /i is spectral with spectrum A. 

Proposition 2.2. [JP98b, lDJ07cj Let fx be a probability measure on R d . A subset A ofR. d 
is a spectrum for fi iff 



Remark 2.3. Note that in equation (I2.2p . we do not have to worry about possible repetitions 
in A. Indeed, if A appears at least twice in (I2.2p . then take x = —A. Since /i is a probability 
measure, /t(0) = 1. So, the sum on the left of (12.21) is at least 2. This would contradict (12. 2p . 
and therefore a A can appear at most once in the sum. 

Affine IFSs. Let R be a d x d expanding integer matrix, i.e., all eigenvalues A satisfy 
|A| > 1. Let B be a finite subset of Z d of cardinality j^B =: N, with G B. We consider 
the iterated function system 



2. Definitions 



(2.1) 




(2.2) 




(2.3) 



r b {x) = R~\x + b), (x G Mr, b G 



B)- 
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Theorem 2.4. |Hut81] There exists a unique compact set Xb C M. d such that 



(2.4) 



X B = \^T h {X B ). 



beB 



Moreover 



(2.5) 



X B 



R~% \b k eB for allkeN 



There exists a unique probability measure ji = jj, R B on R such that 




for all continuous functions f on~R d . 
Hadamard triples. 

Definition 2.5. Let R and B as above. Let L C Z d . We say that (R, B, L) form a Hadamard 
triple if j^L = j^B = N, G L and the matrix 



Assumption. Throughout the paper we will assume that (R, B, L) is a Hadamard triple. 

An easy computation shows that both Definitions 12.11 and 12.51 are closed under taking 
tensor product. For example, the complex Hadamard matrices occurring in (12. 7p include the 
matrices defining the Fourier transform on finite abelian groups, as well as tensor products 
of these matrices. Currently there is no complete classification of all the complex Hadamard 
matrices covered by formula (12.71) . Similarly if a fixed complex Hadamard matrix H is given, 
it is of interest to know all the triples (R, B, L) with the property that H is obtained from 
(R,B,L) via formula (12. 7p . Example 15. II below serves to illustrate these issues. 

The purpose of this paper is to explore ways of building higher dimensional spectral 
pairs from combinatorial algebra applied to lower dimensional systems. Theorems 14.41 and 
Corollary 14.51 are cases in point. We also address the converse problem of factoring higher 
dimensional spectral pairs into products of "smaller" ones. 

We will denote by 5* the transpose of R, S := R T . We will need the following "dual" 
iterated function system 



(2.7) 




is unitary. 





Definition 2.6. Let 




2 



(2.9) 



Then 



(2.10) 



^Wb^x)) = 1. 

leL 
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3. Invariant sets and path measures 



The purpose of the present section is to introduce a setting from symbolic dynamics which 
we will use on a particular given IFS subject to the Hadamard-triple law (Definition 12. 5j) . 
The data of a Hadamard triple, the scaling matrix R, and the two finite sets B, L allow us 
to create a digital filter Wb- We then introduce W^-cycles, and invariant sets. We form a 
compact space Q of infinite words in L, and a W^-path- space measure on Q, see Definition 
13.61 The purpose of this is to allow us to construct candidates for spectra A, and then to 
check the condition in Proposition 12.21 

Invariant sets. For x G M , we call a trajectory of x a set of points 

{o"u,„ ■ • -o-^x \n > 1} 

where {u n } n is a sequence of elements in L such that Wb^c^ ■ ■ ■ a^x) ^ for all n > 1. 
We denote by 0(x) the union of all trajectories of x and the closure O(x) is called the orbit 
of x. If Wb(cix) 7^ for some I G L we say that the transition from a; to 07X is possible. 

A closed subset F C M. d is called invariant if it contains the orbit of all of its points. An 
invariant subset is called minimal if it does not contain any proper invariant subsets. Let 



the space of infinite words over the alphabet L. 

Path measures. For all x G M. d , there exists a unique probability measure P x = 
P X {R, B, L) on Q such that for all /1, . . . , l n G L, n G N 



Q ■= L N = {hl 2 ... \l k e L, for all k e N} 



N 



(3.1) 



P x {{UlU) 2 • ■ • G ft I UJx = h, . . . , UJ n = l n }) 



\\W B {(T k ...aix). 



k=l 



The next three propositions can be found in |DJ07cj . 
Proposition 3.1. For all x G M. d 



00 



(3.2) 



\Kx)\ 2 



Y[W B (S- k x) 



k=l 



For all x G R d and all hl 2 



00 



(3.3) 



P x ({hh...}) 



n^K-- 



a h x). 



k=l 




(3.6) 



WBMaO)M^(s)) = h F (x), (x G M d ), 

2eL 



6 



DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 



and for P x -a.e. oo G 

Proposition 3.3. Let Fi, F2, . . . ,F P be a family of mutually disjoint closed invariant subsets 
ofW d such that there is no closed invariant set F with F H [J k F k = 0. Then 

P x (\jN{F k )\ =1 (xeR d ). 

Recall that the triples (R, B, L) in Definition 12.51 involve two sets B and L. The role they 
play is that they are the beginning of a Fourier duality based on the scaling with the matrix 
R on one side of the duality, and with the transposed matrix S = R T on the other side. The 
discussion below is based on this, and the tools we build come from the iterated function 
system based on the pair (S, L) via formulas (12. 8p . 

Theorem 3.4. [CHR971 Theoreme 2.8] Let K be minimal compact invariant set contained 
in the set of zeros of an entire function h on~R d . 

a) There exists V , a proper subspace of M. d invariant for S (possibly reduced to {0}) ; 
such that K is contained in a finite union 1Z of translates of V . 

b) This union contains the translates of V by the elements of a cycle 

{x , ct^Xq, . . . , o\ m _ x ■ ■ ■ 0^X0} contained in K , and for allx in this cycle, the function 
h is zero on x + V . 

c) Suppose the hypothesis "(H) modulo V" is satisfied, i.e., for all p > the equality 
Cjfci • • ■ OfcpO— er fc / ■ • • o~k> p G V, with k{, fc-el implies h—k^ G V for alii G {1, . . . ,p}. 
Then 

TZ = {x + V, a h x + V, . . . , o im _ x ■■■a h x Q + V}, 

and every possible transition from a point in K Dai ■ • ■ ct^xq + V leads to a point in 
K H ai q+1 ■ ■ ■ a h x + V for all 1 < q < m — 1, where oi m ■ ■ ■ a h x = x . 

d) Since the function Wb is entire, the union TZ is itself invariant. 

A particular example of a minimal compact invariant set is a W^-cycle. In this case, the 
subspace V in Theorem 13.41 can be chosen to be V = {0}: 

Definition 3.5. A cycle of length m for the IFS (cr;)ieL is a set of (distinct) points of the 
form C := {x , a h x , . . . , oi m _ x ■ ■ ■ a h x }, such that a im ■ ■ ■ a h x = x , with l u . . . , l m G L. A 
Wb- cycle is a cycle C such that Wb(x) = 1 for all x G C. 

For a finite sequence l±, . . . , l m G L we will denote by l\ . . .l m the path in fl obtained by 
an infinite repetition of this sequence 

l\ • • • lm • (^1 • • • Imll ■ ■ ■ lm ■ ■ •) 

Definition 3.6. A closed invariant set F is called spectral if there exists a subset A(F) of 
R d such that 

(3.8) P X (N{F))= Yl l/^ + ^l 2 ' (^4 

AeA(F) 

In this case A(F) is called a spectrum for F. 

Proposition 3.7. Suppose (-Fj) i=ljn are invariant sets with the following properties: 
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(i) The set Fi is spectral with spectrum A(Fj); 

(ii) For each 1 < k < n and 1 < i\ < i 2 < ■ ■ ■ < ik < n, the set F ix fl F i2 fl • • • fl F ik is 
spectral, with spectrum A(FjJ fl • • • fl A(F ik ). 

Then 

n 

(3.9) P x ([jN(F l ))= Yl l/^ + A)| 2 

i=l AGU^ =1 A(Fi) 

Proof. First we prove that if i*\ and F 2 are invariant sets then 

(3.10) N(Fi fl F 2 ) = N(Fi) fl N(F 2 ). 

The inclusion "c" is clear. For the converse, let u; G N(Fi) fl N(F 2 ). Then 

rf(^...a Wl (0),F 4 ) -> fori =1,2. 

Suppose, by contradiction, that d(a LOn . . .a^Fi fl F 2 ) does not converge to 0. Then there 
exists a 5 > and a subsequence such that y n := c Uk{n) . . . ct^O and d(y n , Fi fl F 2 ) > 5. Since 
is contained in the attractor Xl of the IFS (ai)i e i, it follows that y n G Xl for all n. Since 
Xl is compact we can find a subsequence such that a^ k{n{p)) . . . a wi converges to some point 
x . But then xq must be both in Fx and F 2 , which gives the contradiction. 

By induction, we can extend (13.101) to any finite number of invariant sets. (Note that the 
intersection of invariant sets is itself invariant) 

We have then 

n n k 

^(U Ar (^)) = E(- 1 ) fe+1 E p *(() N ( F < i )) = 

i=l k=l l<ji<--<ife<n j=l 

n k n 

E(-!) fc+1 E ^(f1^)) = E(- 1 ) fc+1 E E im*+a)i 2 = 

k=l l<h<-<i k <n j=l k=l l<h<-<i k <n AeA(n^ =1 Fi ) 

E^ 1 )^ 1 E E IM^ + A)| 2 = J] |^ + A)| 2 . 

k=i l<ii<-<i fc <n Aen| =1 A(F,;.) Aeu™ =1 A(F) 

□ 

Starting with the iterated function system from the pair (£7, L) in formulas (I2.8p . we have 
built a system of VPs-cycles C. The proposition below gives conditions for being able to 
generate a spectrum A(C) directly from this data. 

Proposition 3.8. Lei C be a WB-cycle. Then C is a spectral invariant set with spectrum 
A(C) = t/ie smallest set A i/iai contains — C t/ie property that SA + L C A. 

Proof. Let C =: {xo, Xi, . . . , x m _i} with x& = o~i k Xk-\ for all G {2, . . . , m} and <7/ m x m _i = 
Xo- With |DJ07cJ Lemma 4.1], we know that N(C) consists of the infinite words in Q that 
end in a repetition of the finite word l± . . . l m , i.e., words of the form u ui . . . uj n l\ . . . l m . 

We proved in [DJ07cj Lemma 4.9] that if u — Uq . . . 0Jkm-ih ■ ■ f° r some k > 0, Ui G L, 
then 

(3.11) P x ({u}) = \ft(x + k c (co))\ 2 , (x G R d )\ 
where 

(3.12) k c (uj) =lu + Suj 1 + --- + S^Ukm-i - S km x . 
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Note that every uj G N{C) is of the form uj . . . uii-m-ih ■ ■ ■ \mi where l\ . . . l m is a circular 
permutation of l\ . . . l m . 

Using f)3.1ip . we have to prove only that the set A(ke) of all numbers kc{oS) with uj of the 
form ujq . . . Ukm-ih ■ ■ -lm, is the set A(C) defined in the statement of the Proposition. 

Since = Xj, we have — = 1% — Sxi for all %. 

First, note that if uj — uj . . . uJkm-ih ■ ■ ■ lm, and uj-i G L, then 

Sk c (uj) + uj-i = w_i + Sujo + ••• + ••• + S km uj km ^ - S km+1 x = 
w_i + • • ■ + S km uj km + S km+1 h - S km+2 Xl = 

_ , , i i okm, , I okm+li i i okm+m—li okm+m^ 

■ ■ ■ — UJ-\ -f- ■ ■ ■ -f- O UJkm ~r >J ti -r • • • -T <J 'm-1 — " ^m-l — 

A; c (cj_i . . .a; fem _iZi . . . / m _i 

This shows that SA(kc) + L C A(kc). On the other hand this calculation shows that any 
kciyj) from the set A(kc) can be obtained (in a unique way!) from a point in the cycle by 
applying operations of the form x i— > So; + I with / G L. This implies that A{kc) = A(C). 

The uniqueness of the operations x Sx + I that lead from a cycle point to a fee (a?) 
comes from the fact that elements in L are incongruent mod S"Z d (see |D J06b[ IDJ07c] for 
more details). □ 

The next result shows how the spectral problem transforms under coordinate changes. 

Definition 3.9. We say that two Hadamard triples (R\, B\, L{) and (R 2 , B2, L2) are conju- 
gate if there exists a matrix M G GLd(Z) (i.e., M is invertible, and M and M _1 have integer 
entries) such that R 2 = Mli.M ; . B 2 = MB X and L 2 = (M T )~ 1 L 1 . 

If the two systems are conjugate then the transition between the IFSs (t^^Bx and (TMb)beB! 
is done by the matrix M; and the transition between the IFSs (<Ti)i€Li an d {o~(M T )- 1 i)ieL 1 is 
done by the matrix (M 7 ")" 1 . 

Proposition 3.10. // (Ri, B 1 , L x ) and (R 2 ,B 2 ,L 2 ) are conjugate through the matrix M, 
then 

(i) r Mbl (Mx) = Mr bl (x), a {MTr i h ((M T )- l x) = (M T )~ 1 (Ti 1 (x) , for all h G B x , h G L x ; 

(ii) W B2 {{M T )- l x) = W Bl (x) for all x G R d ; 

(iii) For the Fourier transform of the corresponding invariant measures, the following 
relation holds: /1b 2 ((M t )~ 1 x) = /1b 1 (x) for all x G M d ; 

(iv) Foruj G Qi = ^et (M t ) _1 cj := ((M t )- 1 uj 1 , (M t )- 1 uj 2 , ...). The associated path 
measures satisfy the following relation: 

pl^m^E) = pi{e). 

(v) If Fi is invariant for the first IFS, then (M T )~ 1 Fi is invariant for the second IFS. 
If in addition Fi has spectrum A(i*i) , then (M T )~ 1 Fi has spectrum (M T )~ 1 A(F 1 ). 



Proof. The results follow by direct computation. 



□ 
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4. Invariant subspaces 

In this section we give two ways of building higher dimensional spectral pairs from com- 
binatorial algebra and spectral theory applied to lower dimensional systems. Theorems 14.41 
and Corollary 14.51 are cases in point. 

Remark 4.1. If V is a rational invariant subspace for S, i.e., it has a basis consisting 
of vectors with rational components then, by [CHR971 Lemme 4.2], there exists a matrix 
M G GL^(Z) as in Proposition 13. 101 that maps V into lR r x {0}, with r = dimV^. Therefore 
we can reduce the study to the case when V = W x {0}. 



Suppose the subspace W x {0} is invariant for S. Then the matrix S has the form 
(4.1) S 



Sx C 
s 2 



(4.2) 



where S\ is r x r, S 2 is (d — r) x (d — r), C is (d — r) x r and is r x (d — r). 
The matrix R has the form: 

A^ 1 



R 



A x 

C* A 2 



and R' 1 



— A~ 1 r'*A~ i A~ 



with A x = S?, A 2 = ST, C* = C T . By induction, 



(4.3) 



Br 







D k A 2 



k-l 



where D k := - ^ A 2 {l+l) C* A~ x 



(k-i) 



1=0 



The set B can be written as 

(4.4) B = {(n, %j ) I i G {1, . . . , Nx], j G {1, . . . , N 2 (t)}} 

with {n, ...,r Nl } = proj Rr . B. 
We have 



X B = {^ J R- fc 6 fe |6 fe G J B}. 



fe=i 



Therefore any element (x, y) in Xb can be written in the following form: 



x 



5>r 



fe=l 



fe=l 



fc=l 



Define 



X i: ={^VnJ^G {!,..., iVj}. 



fe=i 

Let /^i be the invariant measure for the iterated function system 

r n {x) = A 1 - 1 (x + r i ), z6{l,...,JVi}. 
The set X\ is the attractor of this iterated function system. 

Assumption. We will assume that the measure fix has no overlap, i.e., /ii(r ri (Xx) D 

7 r; (.V;)) (I/Or / / ./. 

We will also use the contractions 

a t (x) = S 2 \x + l), {x,leR d - r ). 
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For each sequence u = {ixi 2 . . .) G {1, . . . , Ni} N =: Qi, define x(u) = YlT=i ^i~ fcr V Also, 
because of the non-overlap condition, for /ii-a.e. x G Xi, there is a unique to such that 
x(uj) = x. We define this as uj(x). This establishes an a.e. bijective correspondence between 
Qi and X\, lo <-> x[uS). 

For w = («i«2 ■ ■ ■) E Qi define 

Q 2 {uj) := {Vil,jlVi2,j2 ■ ■ ■ Vin,jn ■■■ 1 h £ { 1 , ■ ■ ■ , N 2 (i fe )}}. 

For u; G fii 

oo 

(4.5) flr(w) := ^D k r ik , and £f(x) := ^(w(ar)). 

k=l 

Also we denote f2 2 (x) := Q 2 (oo(x)). 
For define 



j k G {!,..., iV 2 (ifc)} for all k 



X 2 (x) :=X 2 {u{x)):=lY,A 2 \ k , 3k 

I fe=i 

Note that the attractor Xb has the following form: 

X B = {(x, g{x) + y) | x G X 1; y G X 2 (s)}. 

We will show that the measure /x can also be decomposed as a product between the measure 
Hi and some measures /i 2 on X 2 {oj). 

On fi 2 (u;), consider the product probability measure (jl(u)) which assigns to each r] ik j k equal 
probabilities 1/N 2 {i k ). 

Next we define the measure on X 2 (u). Let r w : f2 2 (u;) — > X 2 (u;), 

oo 
k=l 

Define the measure /J 2 , := fJ>u( x ) := ^{ x ) ° r J(x)- 

In building higher dimensional spectral pairs from lower dimensional systems by factoring 
it turns out that overlap properties in the lower dimensional systems play a critical role. 
This is made precise in the next Proposition. The proof is contained in [DJ07cj : note that 
the more restrictive assumptions used there are not needed for this Proposition. 

Proposition 4.2. Using the notations above, suppose 

(a) The measure fii has no overlap, i.e., fj,i(T ri (Xi) fl r rj (Xi)) = for all i ^ j. 

(b) The numbers N 2 (i) = {rj | (r^, rf) G B} are all the same, N 2 (i) = N 2 , i G {1, . . . , Ni}. 
In particular NiN 2 = N. 

Let a be the shift on a{i\% 2 . . .) = (i 2 i 3 . . .). Let u = {iii 2 . . .) G f^. 

(i) For all measurable sets E in X 2 (uj), 

N 2 (ii) 

2 i=i 

The Fourier transform of the measure satisfies the equation: 
(4-6) p>l(y) = m(S 2 l y,ii)fil {u)) (S 2 l y), 
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where 



, n 2 

m(y, ll ) = —J2e 2m ^ y . 



N 2 . , 



^ii) For all continuous functions on 



fdn= I I f(x,y + g{x))dii 2 x (y)dn 1 {x). 



J yx, y ^ yyx)) a^- 2 ' 

'X B JX 1 JX 2 (x) 

(hi) If Ai is a spectrum for the measure \L\, then 

t>d—r\ 



F(y):= J2 \Kx + Xi,v)\ 2 = [ \fi(v)\ 2 dfi 



x[s) {xER',yE 



A1GA1 

Moreover, if 

Ni 



(4.7) W( y ):=^^T\m(y,i)\\ (y E R d ~ r ) 



then 



(4.8) F(y) = W(S^ l y)F{S^y), and F(y) = J] W(S^ k y), {y E R d ' r 



k=l 

Definition 4.3. We say that p E M. d is an S-period for a function / on K d if 

fix + S n p) = f(x), for all x E R d , n > 0. 

Theorem 4.4. Suppose yo E K d ~ r and the set R r x {y } is invariant. Then W x {0} is 
invariant for S. Using the notations above, there exists (l®, G £ suc/i i/iai o"/o(?/o) — 2/o- 

Assume that the conditions (a) and (b) in Proposition \4-^\ are satisfied. Assume in addition 
that the following conditions are satisfied: 

(i) The measure \L\ is spectral with spectrum A\ . 

(ii) The spectrum Ai has the following properties: 

(4.9) 5(Ax H„}) + LD Ax {-y }. 
and 

(4.10) (Ai, —yo) is an S-period for Wb, (Ai E Ai). 

(hi) (Ax, {rj | i E {1, . . . , A?i}, := {^i I (^i, Z°) e L}) ^ s a Hadamard triple. 

Then W x {y } 2s a spectral set with spectrum 

(4.11) A := |J (L + SL + • • • + S n ~ l L + ^(Aj x {-yo})). 

neN 

Proof. Take x G lR r , then from (I2.10p we deduce that there exists I = (Z?,^) ^ ^ such that 
2/o)) 7^ 0. Then for a;' E W small enough, we have that Wb(ci(x + x',y )) 7^ 0. 
Therefore the transition (x + x',y ) 1— > ai(x + x',y ) is possible. But, since R r x {y } is 
invariant, it follows that W x {y } 3 ai(x + x',yo) = S' _1 (x + x',yo). Subtracting, we get 
S-\x', 0) E W x {0}. Then we can multiply by scalars to see that S'^x, 0) E W x {0} for 
all x E K r , so W x {0} is invariant for S. 
Then we have 

<7,(z, yo) = S" 1 ^ + Z°, y + = (^(x + «?) + Dj( yo + S^(y + l° 2 )) eW x {</ } 
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which implies that cr z o(y ) = Vo- 
Next, we compute for (x, y) G M d : 

W B{°{hil){x,y)) = Yl Yl e 27ri(ri_r * ) - (5 i" 1(a:+ ' l)+I? ? (2/+ ^))+(*^-^',/)-(^ 1 (2/+'g)) 

/l£ii(/§) «l i,i'=lj,j'=l 

= (*) 

But using the Hadamard property in (iii) we have 

J_ _2 7 ri(r i -r<).Sr 1 h _ / 1, * = *' 

So 

1 1 M 

1 2 i jj' 1 »=1 

We obtain the equation 

(4.12) £ ^(^^(^y)) =W r (^(y)) l (?/eR d - r ). 

Using |DJ07cl Lemma 4.1], we have that iV(IR r x {yo}) consists of all the infinite words in 
L N for which the second components are l 2 from some point on. 
We show first that 

oo 

(4.13) P(x,y){{{ui ...cu n ...)\ uj n)2 = l 2 for all n}) = JJ W(aq ■■■a l oy). 

k times 

We compute for all n, by summing over all the possibilities for the first component, and 
using (J33|): 

n 

P{x,y){{{UlU2...)\u) kj2 = l° 2 ,l<k<n})= n W B(<7(l k ,ul 2 ) ■ 4 4 a (h,l,l%)( X 'V)) = (*)■ 

h,i,---,ln,i k=\ 

Using (I4.12p we obtain further 

n-1 

(*) = W(aq---aqy) II ^(^m.^) " ' a (h,i,i 2 )( x ^)) = ■ ■ ■ 

h,l,—,ln-l,l k=l 

n 

= Y[W(aq---aqy). 

k times 

Then, letting n^cowe obtain ( 14.131) . 
We have 

<r${v) = S 2 \y + l° 2 ) = S 2 \y - y ) + S^(y + l° 2 ) = S 2 \y - y ) + y . 
°i°vi°Ay) = S 2 2 (y - y ) + S^iyo + Z°) = S 2 2 (y - y ) + y . 
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By induction 

(4-14) a l o--jaq=S^ k {y-y )+y , (k > 1) 

k times 

Therefore 

oo 

(4.15) P(x,y)({(wi ■■■uj n ...)\ u n>2 = l° 2 for all n}) = J] W{S 2 \y - y ) + y ). 

k=l 

Since M r x {y } is invariant, we have W B (ai(x, yo)) = for all x G W and I = (/i, / 2 ) G L 
with l 2 ^ I 2 °. This implies that 

W(y ) = W^oyo) ^ £ W B ((7 (ll ^ ) (0 lW) )) = X;^(a,(0, l to)) ty S 325 !. 

But then, with (14. 7p . |m(y ,OI = 1 f° r a ^ * e {!)••■ Using the formula for m(y ,i) 

we get that for each i G {1, . . . , Ai} the numbers e 2mVi ' j , j G {1, . . . , A^} are all the same 
(otherwise \m(y ,i)\ < 1). And this shows that \m(y+y ,i) \ = \m(y,i) \ so W(y+y ) = W(y) 
for all y G R d - r . 

Plug this into (14.151) . we get using Proposition I4.2( iii) and the notations there: 



(4.16) P^Mfa ...u n ...)\ u n , 2 = l° 2 for all n}) = J] W{S 2 \y - y )) = F(y - y ) 

k=l 

= + Ai, 2/ — 2/ )| 2 - 

AiSAi 

Next, fix l\, . . . , l n G L. We compute, using (13. 3ft and ( 14.151) 

(4.17) P( x>y )({u \uJi = h,...,uj n = l n ,Uk,2 = 1% for k > n + 1}) = 

W B (a h (x,y)) . ..W B (a ln . . . a h (x, y))Pa ln ...a h (x, y )({u I ^k,2 = 1%, for all k > 1}) = 
WaC^a;,!/)) . ..W B (er/ B . . . a h (x, y)) £ \fl(ai n . ..a h (x,y) + (Xi,-y ))\ 



AieAi 



But 

W B (^(x, y)) y) + (A 1; -y )) | 2 = ^^((z, y) + OM^x, y) + 1 + S(A 1; -y )))| 2 

by pi Wb{s -i {{Xj y) + l + s{Xu _ yo))ms ~\( x , y ) + l + S(Xi, -y )))| 2 

= \jl((x,y) + l + S(\i,-y ))\ 2 . 

By induction, 

W B (a h (x,y)) . ..W B (a in . ..a h (x,y))\jl(ai n . ..a h (x,y) + (Ai,-y ))| 2 = 

\fi((x, y) + l n + Sl n -! + ■■■ + S n -% + S n (X u -y ))\ 2 . 

Using this in (14.171) . we get 

(4.18) P( x>y) ({uj \oj x G L, . . . ,uj n G L,u ki2 = l 2 ,o for k > n + 1}) = 

y) + ln + SL., + ■■■ + S n -% + ETiXx, -y ))| 2 . 

h,...,l n eL,XieA 
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We know by |DJ07cl Lemma 4.1], that N(W x {yo}) consists of all the infinite words in 
L N for which the second component are 1% from some point on. 

By ( 14.91) . the sets L + SL H — • + S n ~ l L + S n (Ai x {—yo}) are increasing with n. Therefore 



P (x , y) (N(W x {y })) 



£ IM^ + A)| 2 . 

\€U n (L+SL+-+S n - 1 L+S n (A 1 x{-y })) 



□ 



In the next result we combine the conditions from the theorem into a form that can be 
used in checking concrete examples. 

Corollary 4.5. Suppose the set W x {0} is invariant. Then it is also invariant for the 
matrix S. Let Li(0) = {li \ {h,0) G L}. Assume condition (a) in Proposition \4 ■ 2\ is satisfied, 
and in addition, the following hold: 

(i) The measure fii is spectral with spectrum Ai. 

(ii) SiAi + Li(0) D Ai, and every X 1 in A 1 is an S\-period for 



Wi(x) :-- 



1 Nl 

-Y 



i=i 



[x e 



Then W x {0} is a spectral invariant set, with spectrum 

A(R r x {0}) := [J (L + SL + • • • + S n ~ l L + S n {A 1 x {0})). 



n>0 



Proof. We use Theorem 14.41 We get that W x {0} is invariant for S, and Z° — in our case. 
We have to check that all conditions are satisfied. Since W x {0} is an invariant set, we have 
that W B (o-i(x, 0)) = for all / G L with l 2 ^ 0. This implies, using ([2TTUD that 

1= W B (a {hfi) (x,0))= W B (S^(x + h),0) = 



N\ N 2 (i)N 2 (i') 



_/y2 



N- 



iieLi(o) 



E E 



i,i'=l i=i i'=i 



E 

ZieLi(o) 



AT- 

E 



i=l 



^ N 2 (i) 27r i r .. ( 5 1 - 1 fa+; 1 )) 



iV 



hGii(O) t,i'=l 
2 

= £ i^^ + zo)! 2 , 

il6Xi(0) 



52i=i ^r 1 ^, and 5 is its Fourier transform. Since 



N Ur i1 



where 5 is the discrete measure 6 
St-^C*) — A, the measure 5 is a probability measure. 

Then, the previous calculation and Proposition 12.21 implies that S'f 1 Li(0) is a spectrum 
for the measure 5. With [DJ07cl Lemma 2.7] we obtain that the numbers iVjs(z) /iV are all 
equal. This shows that condition (b) in Proposition 14.21 is satisfied. Also, since S' ] ~ 1 L 1 (0) is a 
spectrum for 5 this implies that #L%(Q) = N% and (A%, {r 1; . . . , r Nl }, Li(0)) is a Hadamard 
triple. 

Finally, for Ai G Ai we have S n (Ai,0) = (5^Ai,0) so (Ai,0) is an S-period for W B . 
Thus all the conditions of Theorem 14.41 are satisfied and the result follows. □ 
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5. Examples 



The following example is a spectral pair in R 2 whose Hadamard triple is not reducible to 
R x {0} , i.e., to the first coordinate. The Hadamard matrix of the system is of the form 
H ®H where H is the 2x2 unitary matrix of the Fourier transform on Z 2 , the cyclic group 
of order 2. 

Example 5.1. We give an example of an affine IFS that satisfies the Hadamard condition, 
but not the reducibility condition of [DJ07ct Theorem 3.8]. 
Let 



R :-- 



4 
4 



B 













1 




1 





1 


2 


> 


4 




6 



L 








2 




2 










1 





1 


1 


1 


5 



Then the matrix in (12.71) is 



1 



which is unitary so (R, B, L) is a Hadamard pair. 
The subspace R x {0} is invariant. Indeed, we have 



1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


-1 


1 


-1 



W B {x,y) 



_(1 _|_ e ^i2y e 2wi(x+4y) e 2m(x+6y)^ 



-(1 + e 2ni2y )(l + e 27 ™0+ 4 y) x 



And 



a 



(0,0) 



>,0) 



x 



0),cr (2i0) (x, 0) 



,x + 2 



, 0), <T(2,1) (X,0) 



,x + 2 



K 4'4' 



so W^(a (2il) (x, 0)) = H/ B (a (0 , 5 )(x, 0)) = 0, for all x G R. 

However the Hadamard triple is not reducible toRxO because the set L does not satisfy 
condition [DJ07cl Definition 3.1(iii)]: the number of vectors in L that have the second 
component is 2, and there is only one vector that has the second component 1. 

Theorem 5.2. The measure fi = fiB is spectral. 

Proof. We use Theorem I3.4[ and look for candidates for subspaces V such that some finite 
union TZ = U" =1 (xj, yi)+V is invariant and contains some minimal invariant set K. Note that 
since K is minimal, every orbit of every point in K is equal to K. Thus we can take (xi, yi) 
to be in K and also be limit points of some trajectory. So we can assume (x,, yi) G Xl, the 
attractor of the IFS (o";)/ 6 l. 

Using ([23]), we have that X L C [0, §] x [0, §]. 

Since TZ is invariant we must have for all (x, y) G (xj, y^) + V and / G L either o~i(x, y) G TZ 
or Wb(o~i(x, y)) =0. Since X L is not contained in any finite union of translates of a subspace, 
we must have Wb(o~i(x, y)) = for some (x, y) G TZ and some / G L. 

Looking at the formula for Wb we see that WB{x',y') = iff 4y' = 2k + 1 for some 
k G Z or 2x' + 8y' = 2k + 1 for some 6 Z. Thus, the subspaces we are looking for are 
Vi = {(x, y) I y = 0} = R x {0} and V 2 = {(x, y) \ x + Ay = 0}. 

We analyze each subspace separately, and compute what are the possible unions TZ. 



16 DORIN ERVIN DUTKAY AND PALLE E.T. JORGENSEN 

For V\ = R X {0}: since we have Wsi&iixi, yi)) = for some / = (Z 1; l 2 ) this implies that 
4((?/j + Z 2 ) /4) is of the form 2k + 1 for some k G Z. So G Z. And since (xi,yi) G C 
[0, 2/3] x [0, 5/3] this implies that ^ = or ^ = 1. 

For yi = we obtain the invariant set lx{0}. For yi = 1, since we have <7(2,i)(*, 1) = (*, \) 
and Wb{*, 7^ 0, it follows that (*, |) G 7?.. By induction (*, G 7Z, which contradicts 
the fact that 7Z is a finite union of translates of R X {0}. 

Then for V 2 = {(x,y)\x + 4y = 0} we will use a matrix M to conjugate our IFS to another 

[4 _ 1] 

one for which V 2 becomes R x {0}, and we use Proposition 13.101 Take M — n 1 Then 



R = MRM~ 



1 



B = MB = B:- 



L := {M T )- l L 








-2 









-2 










5 


1 




1 














-1 




-5 





5 


2 


5 




6 


5 


20 



1 



We look for 



W s (x,y) = -(1 + e 2 ™(-^)(l + e 2 ^). 

The subspace V 2 = {(4t, -t) \ t G R} is mapped into R x {0} by (M T )~ 1 . 
possible invariant unions 7Z of the form Uj(xj, j/j) + R x {0}. 

We note that C [— |, 0] x [0, y]. Since crj(*, = for some / G Z we get 2- 
for some /c G Z. Therefore must be an even integer in [0, y]. Thus G {0, 2, 4, 6} 

We will use the following notation: if (*, a) is in R and for some I = I2) we have that 
2 • Sy 2 - is not an odd integer, then W B (<Ji(*,a)) 7^ so (* 



Vi+h 



2k+l 



a 



Hyi 

union 71. If y. 



0, then ^ f = 5 







2, then 2A^ 



and this contradicts the fmiteness of the 

1 



... ... ...... jj,.., , , / , 9 -j 2 ") must be in 1Z; we write 

4 • 

s A A A 

4 16 ' ' ' ' 

2 + 2 = lA I A.... H w = 4then4AlAJ_.... „ 

yi = Q then 6 2 1 — > | . . . . In all cases we obtain a contradiction. Therefore V2 does 
not produce invariant sets correlated to minimal invariant sets K as in Theorem 13.41 (one 
needs to apply the conjugation matrix back to (R, B, L)). 

Thus, the only invariant sets we have to worry about are R x {0} and possible W^-cycles. 

First, let us see what the spectrum is for R x {0}. We use Corollary 14.51 We have Ni = 2, 
N 2 = 2, {r 1; r 2 } = {0, 1}. The measure fit is associated to the IFS t (x) = | , T\X = ^4-, so 
it clearly has no overlap. The measure \i\ is a spectral measure with spectrum 



Ai := <^ ^4 fc a fc |a fe G {0,2} 



. k=0 



(This is just a rescale of the first example of a fractal spectral measure given by Jorgensen 
and Pedersen in [JP98bj ). Since A x C Z we see that condition (ii) in Corollary 14.51 is satisfied. 
Hence R x {0} has spectrum 

A(R x {0}) = [j{L + SL + ... S n ~ l L + 5 n (Ai x {0})). 

n>0 



It remains to look for W^-cycles. By |DJ07bl Theorem 4.1] we have to compute the lattice 
r = {7 G R 2 I 7 • b G Z for all b G B}. So if (x, y) G T then 2x, x + Ay, x + 6y G Z so 2y G Z 
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and x G Z, which means that T = Z x |Z. To find the Pv^-cycles we intersect T with 
X L . Since X L C [0, |] x [0, |] it follows that the only candidates for W^-cycle points are 
(0,0), (0,|), (0,1). 

(0,0) is the trivial cycle, but we can discard it because it is contained in R x {0}. It is 
easy to check that the other two are not W^-cycles. 

We conclude that the invariant set R x {0} contains all minimal invariant sets, and there- 
fore, by Proposition 13.31 the spectrum of \i is A(R x {0}). □ 
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